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ABSTRACT 

Charge Transfer Inefficiency (CTI) due to radiation damage above the Earth's atmo- 
sphere creates spurious traiUng in Hubble Space Telescope (HST) images. Radiation 
damage also creates unrelated warm pixels - but these happen to be perfect for mea- 
suring CTI. We model CTI in the Advanced Camera for Surveys (ACS)/Wide Field 
Channel (WFC) and construct a physically motivated correction scheme. This op- 
erates on raw data, rather than secondary science products, by returning individual 
electrons to pixels from which they were unintentionally dragged during readout. We 
apply our correction to images from the HST COSMOS survey, successfully reduc- 
ing the CTI trails by a factor of ^ 30 everywhere in the CCD and at all flux levels. 
We quantify changes in galaxy photometry, astrometry and shape. The remarkable 
97% level of correction is more than sufficient to enable a (forthcoming) reanalysis of 
downstream science products, and the collection of larger surveys. 
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1 INTRODUCTION 

Charge- Coupled Device (CCD) imaging detectors convert 
incident photons into electrons. The electrons are stored 
within a silicon substrate, in a pixellated grid of electrostatic 
potential wells that gradually fill up during an exposure. At 
the end of the exposure, the electrons are shuflfled, row by 
row, to a readout register at the edge of the device. They 
are then counted and converted into a digital signal. 

Above the Earth's atmosphere, a continuous bombard- 
ment of high energy particles makes a harsh environment 
for sensitive electronic equipment. Protons, neutrons, high- 
energy electrons and heavy ions create bulk damage by col- 
liding with and displacing atoms from the silicon lattice. 
The dislodged atoms can come to rest in the interstitial 
space, and the vacancies left behind move about the lat- 
tice until they combine with interstitial i mpurities, such 
as ph osphorous, oxygen or another vacancy l|janesick et al\ 
l200ll h Such defects degrade a CCD's ability to shuffle elec- 
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trons, known as its Charge Transfer Efficiency (CTE; Charge 
Transfer /nefficiency CTI=1-CTE). Electrons can become 
temporarily trapped in the local potential, then released af- 
ter a delay that depends upon the properties of the lattice 
and i mpurities, and the operating tem perature of the de- 
vice (jShocklev fc Readlll952l : lHal]||l952l ). If a few electrons 
are trapped during CCD readout, and held while other elec- 
trons are moved along, they are released as a spurious trail 
(see figure[l} . Regions of the image furthest from the readout 
register are worst affected, because electrons starting there 
encounter the most charge traps during their journey across 
the device. 

Charge Transfer Inefficiency gets worse over time, as 
cumulative radiation damage creates more charge traps. By 
January 2007, when the Hubble Space Telescope (HST) Ad- 
vanced Camera for Surveys (ACS) encountered an electronic 
failure, the instrument had been in operation on-orbit for 
almost 5 years and its Scientific Imagaging Technologies 
(SITe) ST002A CCDs were aheady significantly affected by 
CTI. Similarly, the Wide Field and Planetary Camera 2 
(WFPC2) instrument spent 16 years on-orbit. By the time 
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Figure 1. A typical, raw HST ACS/WFC image, in units of electrons. This 30"xl5" (600 X 300 pixels) region is at the far side of the 
CCD to the readout register, which lies towards the bottom of the page. It was obtained on 15 May 2005, 1171 days after the launch of 
ACS. Upon close inspection, as illustrated in the zoomed inset, the CTI trailing behind (above) objects is manifest. 



it was decommissioned, its truly severe CTI created trailing 
that was readily visible on every exposure. 

Different astronomical measurements are liindered by 
different species of charge traps. Charge traps with a char- 
acteristic release time of a few CCD clock cycles move elec- 
trons by a few pixels, creating short trails that alter the 
apparent position (astrometry) and shape (morphology) of 
faint objects. Charge traps with a release time of many 
clock cycles completely detach electrons from their original 
objects, thus also lowering their measured brightness (pho- 
tometry). Several parametric schemes have been developed 
to estimate and thus correct all these effects at a catalogue 
level. Fitting formulae have been approximated, as a func- 
tion of objects' detector position, date of ob servation and 
flux, for incorrect astrometry with WFP C2 ("Cawl ev et al\ 
120021 '). photometry with AC S (jRiessll2003j') . morp holo gy wit h 
STIS ^Rhodes ef a/.ll2004l ) and ACS (|Rhodes^ ef a^.^^2007^ ■ 
and others. Although such schemes provide first order mean 
correction for a large number of objects, they are inaccurate 
for high precision measurement of individual objects. For 
example, no account is made for the screening of faint ob- 
jects by bright objects slightly closer to the readout register 
(which pre-fill charge traps) or, with extended sources, for 
the dependence upon object size , radial profile, and shape 
(see discussion in lRiessll200Gl and [Rhodes et alll2009l ). 

More precise removal of CTI trailing requires manipu- 
lation of the raw pixel data, using software to move elec- 
trons back to their original locations. Since CCD readout 
is the last stage of data acquisition, this should be the 
first stage of a final d ata reduction pipeline {e.g. CALACS, 
iPavlovskv et al ]|2006l ). Even this will never achieve precision 
at the level of single electrons, since charge trapping and re- 
lease are stochastic, quantum mechanical effects. Ifowever, 
such software has been demonstrated to be much more ac- 
curate than pa rametric solutions for the correction of HST 
STIS imaging (|Bristow l '2003a'; 'Piat ek et al. I l2005l . 2006 & 
2007) and spe ctroscopy (Bristow 2003^), and Chandra ACIS 
spectroscopy (jCrant et o/.l 12004 ). This approach can also 



provide physically motivated model parameters, rather than 
an ad-hoc fitting function to arbitrary parameters. 

In this paper, we develop an empirical but phys- 
ically motivated, pixel-based CTI correction scheme for 
ACS/WFC imaging. We concentrate on traps with release 
times of a few CCD clock cycles, but our algorithm is suffi- 
ciently general to be able to incorporate additional species 
of traps with longer release times if necessary. In we 
build a model of the readout process by measuring the rate 
at which the CCDs' potential wells are filled by electrons, 
and the properties of charge traps within those wells. In 
we describe software to impl ement th i s reado ut model, and 
use the iterative approach of iBristowl (|2003al ) to reverse the 
process. In J?) we apply this to real data and obtain better 
measurements of galaxy photometry, astrometry and ellip- 
ticity in the HST COSMOS survey. In iJS] we discuss possible 
improvements to our algorithm, and summarise our work. 



2 MODELLING THE ACS/WFC CCDS 
2.1 Data and dates 

We primarily base our analysis on the "_raw" data obtained 
during HST cycles 12 and 13 for the COSMOS survey (HST- 
GO-09822, P.I.: N. Scoville). These are 2368 uniform, extra- 
galactic exposures of 5 07 seconds each (|Scoville et aZ.|[2007l : 
iKoekemoer et al\\200^ ). For convenience, we first split each 
exposure into the four quadrants read out through different 
amplifiers, rotating and flipping them to orient the readout 
in the same direction for all. We then multiply the image by 
the calibrated gain (which varies slightly between amplifiers) 
and subtract the corresponding superbias images created by 
STScI. This reverts the images to units of electrons, in ap- 
proximately their configuration on the CCD immediately 
prior to readout (but note uncertainty about the imperfectly 
modeled injection of charge by warm pixels during readout 
in ^2:211 . 
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We have also tested our CTE correction method in a 
less systematic manner on imaging from cycle 14 (HST- 
GO-10496, P.I.: S. Perlmutter). This extends the analysis to 
nearly January 2007, when ACS failed. While this smaller 
survey provided sufficient data to verify an extrapolation of 
our correction parameters, such non-uniform observational 
strategies proved less useful for tracking the gradual CTE 
degradation. 

Our model could also be extrapolated to correct ob- 
servations obtained after servicing mission 4. However, this 
would be a large extrapolation, and the trap density did not 
necessarily continue to rise at the same rate while the in- 
strument was offline. We would therefore recommend new 
calibrations as soon as sufficient images become available. 

2.2 Warm pixels 

A detector's CTI can be measured using First Pixel Re- 
sponse (FPR) tests. Exten ded Pixel Edge Respo nse (EPER) 
tests, or '^'^Fe radiation (jjanesick et al\ l200ll ). FPR and 
EPER tests involve uniformly illuminating a CCD and mea- 
suring deviations from that uniformity in the few pixels near- 
est the readout register and (virtual) pixels obtained by con- 
tinuing to clock charge past the end of the CCD. In addition, 
a radioactive iron source can be used to create a series of 
5-functions on the CCD, all of a fixed energy level. Deficien- 
cies in the observed number of electrons in some events, or 
deviations in their shape, reveal the presence of charge traps 
between the event and the readout electronics. 

A continuous programme of FPR and EPER moni- 
toring has been carried out on the ACS/WFC detectors 
since launch (HST progra ms 9649,10044,10369 and 10732; 
iMutchler fc Siriannill2005l ) . Radiation testing was performed 
prior to launch, but is not possible on-orbit in HST. 

After experimenting on all three types of data ac- 
quired in a laboratory environment, we decided that that 
the parameters needed for our model are most easily 
extracted from radiation tests. Although in-flight ^^Fe 
data are unavailable, cosmic ray damage has also created 
randomly distributed "warm pixe ls" that suffice instead 
ijBiretta fc Kozhurina-Plataisll2005l ). Warm pixels are short 
circuits in the CCD electrostatic potentials used to collect 
charge, and continuously inject spurious charge into the de- 
vice. During a long exposure, warm pixels create single-pixel 
5-functions. Since the rate of charge injection varies between 
warm pixels, these (5-functions have a range of amplitudes. 
In some ways, this is actually more convenient than an ^^Fe 
source. 

To identify warm pixels within each image /, we locate 
2D local maxima, after unsharp masking to eliminate ex- 
tended objects (whether they be resolved galaxies or merely 
a Point Spread Function). One such warm pixel is labeled po 
in figure [T] and contains rie = I{po) electrons. This number 
has been only slightly changed during readout. We measure 
the trail 



(1) 



where p±i are pixels defined relative to the location of the 
warm pixel po in figure[l] and i is an integer from i — 1...9. 
Note that we could have replaced pixels with negative in- 
dices by a locally-determined background level. Our better 
scheme is robust to the accidental inclusion of any extended 



sources (or cosmic rays covering more than one pixel) in our 
catalogue. Although such individual artefacts have a shape, 
they should not have a preferred direction on average, and 
therefore should not bias our measurement of Ti. Indeed, 
we could even have intentionally measured the trailing be- 
hind extended sources. However, in practice, averaging out 
their intrinsic radial profile adds more noise than the extra 
numbers of pixels add signal. 

At the very faint end of our warm pixel catalogue, which 
we push down to rie — 100 electrons so that it only just ex- 
ceeds the background level, our algorithm often finds noise 
peaks instead of true warm pixels. If the noise were due 
to photon shot noise, these pixels would be ideal for our 
purposes, but it is more often read noise - which is added 
after clocking and therefore untrailed. To prevent the dilu- 
tion of our trail measurements, we discard pixels that are 
not flagged as warm in at least half of the exposures. In the 
faintest two bins, this process removes 85% and 50% of the 
catalogue, but quickly becomes negligible for brighter pixels. 
True warm pixels persist through many exposures. 

We avoid "hot" pixels containing more than 76230 elec- 
trons, which is ~ 90% of the full well depth. Saturated pixels 
bleed charge into nearby pixels on the same column, inter- 
fering with the observed shape of the CTI trail. 

Note that warm pixels also inject a low level of charge 
into a column of pixels during readout. This is a relatively 
small complication: the ratio of charge in this line to that in 
the i5-function is equal to the ratio of the clock speed to the 
total exposure time. The data can be corrected along all but 
the worst-affected columns by subtracting superbias frames 
(produced from a zero-second exposure but normal read- 
out). We subtract bias frames before anything else, since we 
are trying to obtain data as it was immediately before read- 
out. However, this is a simpliflcation because (some of) the 
spurious charge was also present during readout, when the 
CTI trails were created. A more complete CCD model might 
also measure the "temperature" of each warm pixel, then in- 
corporate a continuous injection of charge during readout. 
This would be most important for measurements of the den- 
sity of charge trap species with long release times, whose 
extended trail is partly degenerate with a constant line of 
charge injection (as opposed to charge re-emitted from short 
timescale traps for which the original source can be read- 
ily identifled). Even more worryingly, there is (unexplained) 
temporal structure in the injection of charge to some warm 
pixels that is not necessarily reproduced in the superbias 
frames. We therefore decided to incorporate the few percent 
of columns containing the warmest pixels in the CALACS 
bad pixel mask. This is more severe than excluding hot pix- 
els, since the entire rest of the column is masked. 



2.3 Model assumptions and parameters 

To emulate the process of CCD readout, we first need to 
model the collection of electrons within a pixel's poten- 
tial well. In theory, the expected locations of a given num- 
ber of electrons can be calculated from a model of the 
CCD, by solving P oisson's equation within the approp riate 



CUU, by solvmg F oisson s equation withm the approp riate 
applied potent ial. Hardv. Murowinski fc DeenI l| 19981 ) and 
ISeabroke et al\ ()2008h find that, for tvpical CCDs, electrons 
are contained within a well-defined volume - outside which 
their density falls rapidly to zero. Adopting this result, we 
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shall regard a given number of electrons as filling a three 
dimensional pixel up to a specific height. 

We then need to model the charge traps. We assume 
tha t none are fi l led du ring inte gration. Howev er, following 
e.o. lGrant et al\ |200i) (but not iBristowll2003al '). we simplify 
the problem of charge capture as soon as readout begins by 
Eissuming that unoccupied traps instantly capture any elec- 
trons in their vicinity {i.e. all traps inside the well-defined 
volume get filled). In reality, there is an exponential decay 
time needed for traps to capture an electron, but this is typ- 
ically much shorter than the electron dwell time with in one 
CCD clock cycle ()Hardv. Murowinski fc DeenI ligOSl ). The 
rapid and efficient capture of electrons, even in locations 
where the electron density is low, further delineates the 
quantifiable volume of a pixel occupied by a given num- 
ber of electrons. It also means that any charge traps physi- 
cally located near the bottom of a pixel (below the ambient 
background level) reach the steady state of being perma- 
nently occupied, even during readout. In our model such 
traps never affect an image, so we can never probe them; 
but then, we never need to. 

Each filled trap releases its electron after a time de- 
termined by r, the e- folding time of an exponential decay 
IShocklev fc Readlll952l : lHalllll952l ) . If the trap is inside the 
volume currently occupied by electrons, it will immediately 
capture a new one; if it is above, the electron will be pulled 
down and added back to the cloud (which thus expands and 
may even be exposed to new charge traps). We concentrate 
on species of traps with capture times r of a few CCD cloc k 
cycles (3212/xs in the parallel direction: ISirianni et alkoo4 ). 
Traps with long release times affect the leading edge of the 
first bright object in a CCD column, but then reach a full 
steady state (at least while the rest of the object passes by) 
and only mi nimally affect its ph otometry, astrometry and 
morphology (|Rhodes et al\ \2 00^ . For a fixed clock speed, 
r can be conveniently expressed in units of the number of 
pixels that have passed since capture. 

If the charge traps are evenly distributed, we need mea- 
sure only the mean density of each specie^. In ti2.5l we 
indeed find the distribution of traps to be uniform in the 
y direction; we have no reason to suggest that the distri- 
bution is not similarly uniform in the other directions. In 
addition, our approach will internally parametrize out any 
smooth density gradients in the vertical direction. 

We assume that each charge trap always captures a sin- 
gle electron. Real traps might be able to hold several, or oc- 
casionally capture none. If they hold several, and especially 
if the trap densities are low, simultaneous release of many 
electrons could produce correlated spikes in the trails. Trap- 
ping mechanisms are not yet sufficiently understood to know 
whether electrons are released singly or in groups (C. Bebek, 
priv. comm.). We therefore simply assume that any multiple 
(or partial) occupancy of traps is equivalently modelled by 
an the effective density of single-electron traps. 

Finally, we ignore CCDs' three-phase clocking cycle, by 
which a cloud of electrons is transferred from one pixel to 



^ This may not be acceptable for radiation-hard ened detectors 
in fut ure spacecraft that suffer very few traps (see lDawson et all 
l2008h . In that regime, techniques like pocket pumping would be 
ideally suited to measuring the locations of individual traps. 
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Figure 2. The mean CTI trail behind 63901 warm pixels con- 
taining between 3234 and 76230 electrons each, from 20 exposures 
obtained 1171 days after the launch of ACS. Error bars show the 
Icr scatter from measurements behind the whole population of 
warm pixels. The dashed curve shows the best-fit double expo- 
nential model. 



the next (jjanesick et aZ.|[200lh . The transition period ex- 
poses the electron cloud to additional regions of silicon, and 
additional charge traps. However, with our assumption of 
instantaneous charge capture, extra charge traps in regions 
between pixels can again be modeled to first order as an 
increase in the effective trap density. These traps are a frac- 
tion of a pixel nearer to (or further from) the readout reg- 
ister, but horizontal translations of the exponential release 
curve are also degenerate with a change in normalisation. 
The only remaining finesse is that the width (and therefore 
the height) of the electron cloud changes during the three- 
phase clocking cycle, making the cloud dwell temporarily at 
different heights within the silicon and therefore being ex- 
posed to a different number of charge traps. We ignore this 
second-order effect, but incorporating a model of the full, 
three-phase clock cycle might be a useful exercise for future 
work. 

Our model features the following free parameters: 

• number of (relevant) species of charge trap, 

• characteristic release time of each species of charge trap, 

• density of each species of charge trap 

(or a functional form for their distribution if they are not 
uniformly located throughout the CCD), plus the 

• height occupied by a cloud of electrons in a pixel's po- 
tential well {i.e. the number of charge traps they are exposed 
to), as a function of the number of electrons. 

As we shall now demonstrate, warm pixels enable us to mea- 
sure all of these quantities from on-orbit science exposures. 
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2.4 The characteristic release times of different 
species of charge trap 

To obtain maximum signal to noise, we begin by studying 
tiie briglitest warm pixels, furthest from the readout regis- 
ter and in data obtained near the chronological end of the 
COSMOS survey. We select 20 images taken on 15 May 2005, 
1171 days after the launch of ACS, and measure the mean 
trail Ti around all warm pixels between 1634 and 2039 pixels 
from the readout register, and containing between 3234 and 
76230 electron^. The mean trail is shown in figure [21 

Using Levenberg-Marquardt least-squares minimiza- 
tion, we find that the mean trail is well-fit by a sum of two 
exponentials. 



Aie 



-i/Ti 



+ A2e 



(2) 



where Ai = 327±6 e', n = 1Q.4±3.2 pixels, A2 = 108±3 e" 
and T2 = 0.88 ±0.2 pixels. At the — 83°C operating temper- 
ature of ACS/WFC, these correspond precisely w i th the re- 
lease times of trap species found by iHopkinso 3 ll200lli m 
a Marconi CCD47-20 n-channel detector at 0.34 eV and 
0.31 eV below the conduction band. The origin of these par- 
ticular traps is unknown, but they have been speculated 
to be E-centres (Phosphoro us- Vacancy compl exes) with im- 

J iurities of either c arbon l|Hopkinsonl 119911 ) or hydrogen 
Tokuda fc l"toll2000l ). Pure E-centre charge tra ps at . 44 eV 
exist in both the Marconi and ACS/WFC (|Jonesl l200(il ) 
CCDs, with a characteristic release time of te ^ 300 pix- 
els in the Marconi. All ot her charge traps in the Marconi 
CCD have r < 0.01 pixels (|Hopkinsonll200ll ). so this 1-to-l 
correspondence is unlikely to be coincidence. 

However, allowing a third exponential prevents our fit 
from converging, with positive A3 but T3 iterating towards 
infinity. Fitti ng slowly de caying exponentials is a standard 
problem (e.g. ICoveilboOSD . because they are nearly degen- 
erate with a constant. A line of charge above a warm pixel 
could represent either a species of trap with a very long re- 
lease time or contamination from continued charge injection 
during readout (c./. We therefore restrict our analy- 

sis to the two species of traps with shorter r, which also 
have a much more profound efltec t on object astrometry and 
morphology ([Rhodes et fflZ.|[2009i ). Henceforth, we shall fix 
the release times of the charge trap species at ri = 10.4 
and T2 = 0.88 pixels, and also fix their relative densities to 
be 3.0;][f|. Fitting to earlier ACS data yields consistent (but 
noisier) values of r and the relative density. 



2.5 Height of pixel occupied by electrons and the 
density of charge traps 

The total number of electrons that have been displaced from 
a packet originally containing rie electrons is the integral 
under the curve 



1 



A2 



5I/T2 



(3) 



Since we are assuming that any traps surrounded by free 
electrons are instantly filled, this number of trapped elec- 
trons equals the (efltective) number of charge traps Uq in the 
volume Y traversed by electrons between their original pixel 
and the readout register. 

The volume V is the product of the distance y to the 
readout register, the width of a pixel, and the height within 
a pixeI3 filled by rie electrons, above a background level h. 
Warm pixels are available with a range of y and and 
in images with a range of backgrounds h (this varies around 
51±9 e~ due to the an gle between the telescope and the sun: 
iLeauthaud et a/.l2007l ). By counting the number of electrons 
trailed behind a variety of warm pixels, we directly measure 
nq{y, Tie, b). Note that the interpretation of electrons filling 
pixels to a given height is instructive but not necessary: this 
function is precisely what we require. 

Figure [3] shows the mean CTI trail behind warm pixels 
at varying distances from the readout register and of varying 
warmth. The top-right panel reproduces figure [5] The other 
panels show trails behind different pixels within the same set 
of 20 images. Although a busy plot, it has been simplified 
for illustration. For the real analysis, the number of bins was 
doubled in both y and rie directions, and an additional di- 
mension, with five bins of varying background level, has been 
marginalised over here. Towards the left hand side, the trails 
may be systematically underestimated due to the inclusion 
in our catalogue of read noise peaks. However, the efforts in 
i|2.2l to exclude them have kept this small, as demonstrated 
by the still-recovered shape of the trails. Each trail is fitted 
with double exponential decays ([5}, the integral under which 
gives the total number of trapped electrons, and hence the 
number of exposed traps Uq. 

Figure|3]shows the total number of charge traps exposed 
to each set of warm pixels for which a panel was allocated in 
figureO Solid lines join the data points measured at constant 
distances from the readout register (rows in figure [5J. These 
data are well modelled by 



= Pq 



maxjwe — d, 0} 



max{fe — d, 0} 



(4) 



where a, (5, d and pq are to be measured. Parameter w is 
the full well depth, which we assume to be fixed at 84700 e~ 
even th ough it varies b y a few thousand electrons across the 
device (|GiUilandll200i ). and d is the depth of the supple- 
mentary buried channel or "notch" in the CCDs, which we 
fit. The notch is a small region at the bottom of a potential 
well used to gather electrons when only a small number are 
present, much like a narrow channel cut into the bed of a 
water drainage canal. In our model, the notch has zero vol- 
ume, so it merely adjusts the starting point of the power- law 
increase in height governed by a. 

We again use Levenberg-Marquardt least-squares min- 
imization to fit the free parameters. We find that /3 = 
1.01 ± 0.01, consistent with a uniformly increasing number 
of traps with increasing distance from the readout regis- 
ter. Since we would expect this behavior, we explicitly fix 



^ These values are chosen so that more bins can be added later, 
equally spaced in distance from the readout register and log(ne). 
^ There is no reason why we would expect this to be an integer. 



* This assumes that the height of an electron cloud changes by 
only a small amoun t during readout; for further discussion, see 



[Rhodes ef"all |[2009l') . 
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Figure 3. CTI trails behind a range of warm pixels, containing different amounts of flux (increasing left to right) and at greater distances 
from the readout register (increasing bottom to top). Figure |2] is reproduced exactly in the top right panel. The smooth dotted lines 
show fits of double exponential function |2] assuming a fixed ratio of the densities of different species of charge trap, and hence the same 
shaped curve in each case. Open circles depict negative values. These measurements can be used to probe the distribution of charge 
traps on the CCD, and the rate at which a pixel's potential well fills up as a function of the incident flux, thereby exposing electrons to 
additional traps. The mean distances from the readout register for each row, and the flux levels for each column, are provided in figured 
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Figure 4. The total number of charge traps nq exposed between 
a warm pixel and the readout register, as a function of the number 
of electrons rie in that pixel. If traps are uniformly distributed, 
this is proportional to the height within a pixel filled by rie elec- 
trons. Error bars are underestimated, because they do not include 
contributions propagated from uncertainty in r. The solid lines 
join data from warm pixels 212, 618, 1024, 1430 and 1836 ± 203 
pixels from the readout register (bottom to top). The cumulative 
number of traps increases linearly with that distance, implying 
that the traps are uniformly distributed throughout the CCD. 
The dashed curves show a model in which the electrons occupy 
zero volume within the CCD notch (the depth of which is a free 
parameter), and increase as a power law of the number of elec- 
trons above it. 



/3 = 1 and then obtain best-fit values of a = 0.576 ± 0.013, 
d = 96.5 ± 2.0 e" and = 0.544 ± 0.008 pixer\ split be- 
tween the two species of charge trap in a ratio of 3.0:1. Note 
that these errors do not include the uncertainty propagated 
from the measurement of the charge trap release times and 
density ratio. 
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Figure 5. The accumulation of charge traps in the ACS/WFC 
CCDs between 9 September 2003 and the end of 2005. The data 
points and left hand axis show measurements of the total den- 
sity of charge traps, split between the two species in a ratio of 
approximately 3.0:1. Horizontal error bars show the total time 
during which it was necessary to average data to obtain sufficient 
signal to noise. The histogram and right hand axis show the num- 
ber of exposures taken as part of the COSMOS survey, which we 
analysed to obtain these figures. The dotted vertical lines show 
the times when the ACS CCDs were annealed (these are also 
listed in table lAll l. 



density on launch of p° = 0.037 ± 0.001 traps per pixel (al- 
though the formal error on the latter is probably an under- 
estimate, since the extrapolated value is very sensitive to er- 
rors in the gradient). We also tried fitting a sawtooth model 
in which a fraction of the traps were r emoved during their 
first subsequent anneal (as described bv lCox fc CottinghamI 
(2001), the temperature is raised to 13° C about once a 
month; the exact dates durin g cyc les 12 and 13 are reco rded 
in the appendix). iRiessI (|2002l ) and lSirianni et al\ (|2004l ) con- 
firmed that this successfully reduces the density of hot pix- 
els. However, the fraction of removed charge traps iterated 
to ze ro during our fit: recoverin g a linear model (consistent 
with iMutchler fc Siriannillioosi ) and suggesting that the an- 
neals did not significantly improve CTE. Indeed, while hot 
pixels and charge traps are likely both a product of the same 
non-i onisin g energy loss and should accumulate at the same 
rate, iRobbinsi (^2000 ) find that traps are only effectively re- 
moved by annealing at temperatures above 150° C for E- 
centre traps and 330° C for divacancy traps. 



2.6 Growth rate of charge traps 

ACS was installed on 7 March 2002, and cumulative radia- 
tion damage since then has gradually created more charge 
traps. We can use the COSMOS data, which was acquired 
uniformly over a long period, to track CTI degradation. 
Assuming that charge trap release times and the relative 
densities of the two species remain unchanged, we use the 
warmest pixels furthest from the readout register (as in tj2.4p 
to measure the density of charge traps pq at twelve addi- 
tional times during the survey. 

Figure [S] shows the measured increase in charge trap 
density over time. This is fit by a constant accumulation of 
(4.34 ± 0.13) X 10~* traps per pixel per day, and an initial 



3 CCD READOUT SOFTWARE 
3.1 Forward operation 

We now describe a software imitation of the hardware CCD 
readout, using the model developed in Sj2l Just like the 
real version on HST, our software uses a fairly simple al- 
gorithm to shuffie charge to an adjacent pixel, and this is 
repeated many times to complete a full readout. For con- 
venience, our code has a user interface written in IDL, but 
for speed, the core is written in Java. Even so, the software 
algorithm is much slower than an on-board hardware read- 
out, taking 25 minutes per 4096 x 4096 ACS/WFC image on 
a single 2 GHz process or, compared to 2.5 minutes on HST 
dPavlovskv et ttZ.]|2006l ). However, it is trivially parallelisable 
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for multiple exposures, and need ever be run only once on 
any given exposure. 

Since the mean density of charge traps pq is known, but 
their exact locations are not, traps are scattered at random 
3D locations within the model silicon wafer. Importantly, 
each trap is allocated a particular vertical height in some 
detector pixel. The occupancy of each will be continuously 
monitored. Electrons are then added into the pixel array, in 
the numbers we think they were at the end of an exposure. 

A cloud of We electrons is assumed to fill a pixel to a frac- 
tional height (max{ne — d,0} /w)°' . Any charge traps inside 
the occupied volume immediately capture one electron. All 
free electrons are then transferred instantaneously to their 
adjacent pixel (or to the serial register for the last pixel in 
parallel array, or to the preamplifier from the last pixel of 
the serial register). As discussed in ^ we do not attempt 
to model the three-phase clock cycle of a real CCD. Finally, 
electrons inside traps are allowed to decay with probability 
l — e~^^^. Released electrons are returned to the free electron 
pool, and the process is repeated. 

Shot noise in this model created ragged trails behind 
faint objects and amplified background noise fluctuations. 
We introduce two numerical schemes to ameliorate this. 
Firstly, to counter the random release or non-release of elec- 
trons from charge traps, our model traps are allowed to con- 
tain fractional numbers of electrons. In this scheme, each 
electron is gradually released in a floating-point exponen- 
tial trail, until the trap is refilled to a whole electron by the 
passing of a new bright pixel (removing only a fractional 
number of electrons from the free pool). After readout, the 
number of electrons in each pixel is rounded back to an 
integer value for storage. Secondly, to counter the necessar- 
ily random locations of model traps, we can run the code 
multiple times with different random seeds and average the 
results. To save CPU time, this is equivalently implemented 
by increasing the density of traps by a factor of 3 (or any 
other number specified as an optional input parameter), each 
of which is only allowed to contain a third of an electron. 
This most improves the sampling of traps in the thin slice of 
the CCD above the image background level that profoundly 
affect faint sources. With both schemes, the change in back- 
ground noise level during readout becomes less than 1%. 

3.2 Reverse operation 

Moving electrons back to where they belong, and thus re- 
moving the CTI trails, requires the image mapping to be in- 
verted. This is not analytically possible. However, the trail- 
ing represents a small perturbation around the true image, 
so it can be achieved usin g the forward algorithm, vi a an iter- 
ative ap proac h. Following Bristow fc A lcxov (200^). lBristowl 
l|2003al ') and iBristow. Kerber fc Ros~al (|2005), table [1] de- 
scribes a way to obtain an image that, after a software read- 
out, reproduces the data actually downloaded from HST. 
This is the desired, corrected image. Images (A), (C), (E), 
etc can all be obtained and, after several iterations, converge 
to that ideal. 

In practice, excellent results are obtained after only a 
single iteration. Indeed, since the CCD filling and readout 
models are unlikely to match the real hardware better than 
0{S^), only the first iteration is even useful. To speed the 
correction of the 2368 COSMOS exposures, we therefore re- 



True image 


/ 


Not available 


Downloaded from HST 


1 + 


5 


(A) 


After one extra readout 


1 + 


25-1-52 


(B) 


(A)+(A)-(B) 


I ~ 


52 


(C) 


After another readout 




(D) 


(A)+(C)-(D) 


1 + 


53 


(E) 


After another readout 


1 + 


5 + 5^+ 54 


(F) 


(A)+(E)-(F) 


I - 


54 


etc 



Table 1. Iterative method to remove CTI trailing, using only a 
forward algorithm that adds trailing. The true image / is desired 
but not available. Only a version with (a small amount of) trail- 
ing, I +5, can be obtained from the telescope. However, by running 
that image through a software version of the readout process, and 
subtracting the difference, deviations from the true image can be 
reduced to 0{6^). Successive iterations further reduce the trails 
until an image is produced that, when "read out" , reproduces 
data arbitrarily close to those obtained from the telescope. This 
is the corrected image. 



stricted our correction of the larger data set to a single itera- 
tion. Figure|5]shows the image from figure[T]after correction, 
using one iteration. 



3.3 Additional features of the software 

Our CCD readout software has several additional fea- 
tures that we did not find necessary for CTI correction in 
ACS/WFC, but which may be useful for other instruments. 
We document these features here for completeness and fu- 
ture reference. 

Although we used only two species of charge trap, the 
code is capable of handling many (currently, up to four) 
different species, each with a density pq and characteristic 
release time r. With additional measurements of flux loss 
due to long-lifetime charge traps, this could be used to more 
accurately correct photometric measurements of astronomi- 
cal sources. 

The code can shuffle electrons both in the parallel direc- 
tion and along a serial register. In this mode, the execution 
time is doubled. Charge traps in the serial register, with 
characteristic release times similar to the serial clock speed, 
can create trails at 90° to those desc r ibed by [T] However, 
it is well established {e.g. lRie"s3 l2003l : iMutchler fc Siriannil 
l2005h that the serial CTE in the ACS/WFC CCDs is very 
high, and that trailing is almost exclusively in the paral- 
lel direction. Indeed, we attempted to measure serial CTI 
trails but found a signal consistent with zero. We therefore 
neglected this entire mechanism. 

By default, the code places model charge traps ran- 
domly within the model CCD. These are not at the same 
locations as the charge traps in the real CCD, which intro- 
duces additional noise during the correction. This is amelio- 
rated by increasing the density of traps but decreasing their 
capacity as described in >i3.1l Howeve r, radiation-har dened 
detectors in future cameras {e.g. Da wson et a/.ll2008l ) may 
contain sufficiently few traps for this noise to become sig- 
nificant. If their individual locations could be measured on 
orbit, for example via pocket pumping, the code can place 
model traps at the correct locations. 
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Figure 6. Top: The HST ACS/WFC image from figure[T] after CTI correction, in units of electrons. Bottom: Difference image. 



4 TESTING THE CORRECTION 
4.1 Warm pixel trails 

Figure [7] shows the residual CTI trails measured from cor- 
rected versions of the images used for figure [3] Note the 
reduction in the level of all trails by factor of about 30. This 
reduction is matched or bettered during all four of the times 
(the first, sixth, seventh and thirteenth data points in fig- 
ure [Sjl when many COSMOS exposures happened to be ob- 
tained during a single day, and a more precise measurement 
could be made. This demonstrates that we have achieved an 
approximately (1 — 1/30) = 97% precision in our trapping 
and readout model, at all locations on the CCD array and 
for traps at all heights within the CCD. 

If the images are instead corrected using two iterations 
of the method in table [TJ the data points in figure [7] move 
only within their error bars, justifying our use of only one 
iteration for the sake of speed. 



4.2 Galaxy photometry measurement 

We now stack the corrected images to create 577 "_sci" 
science frames, following the pipeline of iKoekemoer et all 



2007^■ which makes use of the Drizzle/Multidriz zle software 



Fruchter fc Hookll2002l : IKoekemoer et al\ \200'^). There are 



four exposures per pointing, dithered almost along a straight 
line (two exposures separated by 6.1" at an angle of 91.2° 
from the coordinate frame and another pair, separated by 
3.1" from the first, at an angle of 85.4°). These are trans- 
formed into a projection of the pixel array onto the sky, with 
oversampled 0.03" pixels that are aligned with the instru- 
ment coordinate frame to keep the two readout registers at 
the top and bottom. The dithers complicate the issue of how 
far each galaxy is from the readout register, because it was 
in different areas of the focal plane in different exposures. To 
roughly divide the catalogue into objects at final positions 
(a;, y) that have traversed different numbers of charge traps, 
we calculate 

^transfers = max {0, 2048 - |6046 - 1.67j/ - 0.0614j;|} (5) 

and separately consider objects with values of ntiansfcrs in 
the ranges 0-475, 475-950, 950-1425, 1435-1900 and 1900- 
2048. The final bin also includes most of the area of the 
oversampled images in which only data was available from 
only three of the four exposures, due to the gap between the 
CCDs. 




Figure 7. Measurement of the CTI trails from the same images used in figure \3\ after correction. The y axes are logarithmic; open 
circles and dotted (rather than dashed) lines show negative values. 
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Figure 8. Changes in galaxy piiotomotry during CTI correc- 
tion, as a function of the measured magnitude and detection 
S/N in the corrected images. A fractional change is calculated 
as (old— new)/new. The five lines connect subsets of galaxies be- 
tween approximately 0-475, 475-950, 950-1425, 1435-1900 and 
1900—2048 transfers from the readout register; the photometry is 
measured on stacked exposures and the exact number of transfers 
varies slightly between dithers. The points show the median value 
in the top panel, since the distribution of points is highly skewed, 
and the mean value in the bottom panel, which is less so. Error 
bars show 68% confidence limits on the mean. 



Figure [8] shows the changes to galaxy photometry dur- 
ing CTE correction, as meas ured on original versus c orrected 
images by So urce Extracto r l|Bertin fc Arnout3ll996l V As ex- 
pected from iRiessI (|2003l ). we find that the most affected 
sources are faint galaxies far from the readout register. The 
amoun t of change is i n hne with the extrapolated predic- 
tion ofH iess fc Mackl (|2004l '). Note however that this should 
not be interpreted as the error in point-source photometry 
on uncorrected images. In particular, our correction scheme 
does not currently include any species of charge traps with 
very long release times, which may steal additional flux from 
sources. Our resolved sources also contain additional elec- 
trons in extended wings that, according to our well filling 
model, lie near the bottom of the pixel potential, occupy a 
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Figure 9. Change in galaxy astrometry during CTI correction, 
in the direction of the readout registers. Readout registers are 
located at opposite edges of the array, with pixels from each half 
being read out at their nearest register. Solid lines connect sub- 
samples of galaxies between magnitude limits of F81AW = 17—19, 
19-21, 21-23, 23-25 and 25-27. For the sake of clarity in this fig- 
ure, the points show medians in each bin. The error bars show 68% 
confidence limits, and the mean values behave consistently around 
these bounds, but noise renders the overall behaviour less clear. 
Note that faint galaxies are also typically smaller than bright 
galaxies, unlike the point sources used in previous studies. 



relatively large volume of silicon and are exposed to a greater 
number of charge traps. 

Photometric measurements can also be expressed in 
terms of the galaxies' detection signal to noise. The inter- 
pretation for galaxies with S/N< 20 is complicated by two 
issues. Firstly, a selection bias in which sources whose de- 
tection signal to noise increased greatly were not included 
in the pre-correction catalogue. Secondly, the correction is 
a form of unsharp masking, and noise in adjacent pixels 
becomes less correla ted during correction. As discussed in 
[Rhodes et oil l|2009l ). this can actually lower the detection 
S/N of faint sources. 



4.3 Galaxy astrometry measurement 

The change in galaxies' positions during CTI correction is 
shown in figure U The spurious shift is 2.42 ± 0.17 x 10"^" 
per (0.05") pixel transfer for galaxies between magnitudes 
17 and 19, and 0.57 ± 0.24" per transfer for 25*'' to 27*'' 
magnitude galaxies. 

The distance depends as expected upon distance from 
the readout register, but is initially su rprising as a fu nction 
of galaxy flux. Using very similar code. lBristowl l|2004 ) found 
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Figure 10. Change in galaxy FWHM size during CTI correction. 
Solid lines connect the same samples as in figure |8] Points show 
mean values within bins and 68% confidence limits. 



bright point sources to be less affected than faint ones. Our 
own code also moves a given object less if its flux is increased. 
However, bright resolved sources also tend to be larger than 
faint ones. Astrome try is most affected by charge traps with 
short release times (jRhodes et alSwO^ ). which are now able 
to capture, move and release an electron several times within 
a large galaxy. Furthermore, the non-linear well filling model 
ensures that even electrons in the low-level, extended wings 
are exposed to most of the charge traps seen by the core. 
Multiple trapping does not remove flux from a galaxy, so 
should not impact photometric measurements. However, this 
process breaks the standard argument whereby charge traps 
saturate, with empty traps being newly exposed only by rel- 
atively large increases in flux. Subtle efects like this highlight 
the power of a pixel-based method: they are unlikely to be 
matched in a parametric scheme acting at a catalogue level. 



4.4 Galaxy shape measurement 

Figure [10] shows the change in galaxies' full width at half 
maximum size during C TI correction, as measured by Source 
Extractor (|Bertin fc Arnouta.1996. '). The CTI trails had spu- 
riously enlarged objects, which are shrunk slightly when the 
flux in the trails is pushed back into the core. The effect in- 
creases as expected with the number of charge transfers, and 
is most pronounced in faint galaxies, which are intrinsically 
smallest, and also encounter the largest number of charge 
traps per electron during readout. 

Galaxy ellipticities are a particularly interesting case 
because the true value of ellipticity is known, at least sta- 
tistically. If there is no preferred direction in the Universe, 
galaxies cannot preferentially align in any direction and their 
mean ellipticity ought to be zero. However, the CTI trails 
coherently elongate galaxies in the readout direction. 

In practice, high precision measurements of galaxy 
shapes are made difficult because they are imprinted 
with the shape of the telescope's PSF. Furthermore, the 
ACS/WFC PSF changes during orbit, as HST expands or 
contracts due to thermal shifts. Most inconveniently, the 
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Figure 11. Component of mean galaxy ellipticities in the direc- 
tion of the readout registers (e cos 28, where e is the ellipticity 
and 8 is the angle between the major axis and the line joining op- 
posite sides of the CCD) before and after CTI correction. Values 
are shown after correction for convolution with the PSF, whose 
shape is imprinted on the galaxies', but this correction increases 
the scatter. Solid lines connect subsamples of galaxies within mag- 
nitude limits of F8UW = 22-23, 23-24, 24-25, 25-26 and 26-27. 
Points show mean values in bins and 68% confidence limits. 



predominant PSF patterns happen to align themselves with 
the readout dir ection and also ten d to align galaxies in 
that direction (|Rhodes et al\ [20o3). Fortunately, a great 
deal of software has been developed for effectively "de- 
convolving" galaxy ellipticities from the PSF, in order to 
me asure weak gravitational lensing. W e use the algorithm 
by Rhodes. Re frcgicr & Groth (2000), as impleme nted by 
iLeauthaud et al. (2007, ) and Massey et al\ (|2007al ). to re- 
cover the underlying galaxy ellipticities. 

The top panel of figure [TT] shows how galaxies appear 
spuriously elongated in the direction of the readout regis- 
ters before CTE correction. The faintest galaxies are worst 
affected, although the elon gation is also s een in the bright- 
est galaxies. As shown in [Rhodes et all (|2009l ). the effect 
is also worst for galaxies intrinsically elongated perpenicu- 
lar to the readout direction. Neit her such effect is mo delled 
by the parametric correction of [Rhodes et al\ (|2007[ ). The 
mean spurious ellipticity of galaxies in the faintest two bins 
increases as approximately (24.3 ± 2.9) x 10~® per (0.05" 
pixel) transfer and (7.7 ±2.0) x 10~® per transfer. The noise 
is predominantly due to the process of PSF deconvolution. 
After CTI correction, as shown in the bottom panel, this ef- 
fect has been successfully reduced to (0.2 ±3.1) x 10~® and 
(2.2±2.1)xl0"^ The mean shear is similarly reduced by an 
order of magnitude from (1.5 ± 0.1)% to (0.2 ± 0.2)% in the 
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faintest bin, and from (7.4±1.1) x 10"'' to (0.9±1.1) x 10"'^ 
overall. 



5 DISCUSSION &: CONCLUSIONS 

We have corrected Charge Transfer Inefficiency trailing in 
ACS/WFC images by more than an order of magnitude. As 
demonstrated by the improvement of HST COSMOS survey 
images between figures |3] and [T] and in figure 1111 we have 
achieved a 90-99% precision in our model, at all locations 
on the CCD array and for traps at all heights within the 
CCD. Our method works at the pixel level, moving individ- 
ual electrons back to locations from which they were dragged 
during readout, and model parameters correspond to real, 
physical properties of the device. This is far bet ter than 
the ad -hoc, parametric correction described bv Rho des et al\ 
(|2007l ). Since the trails are created during the final step of 
data acquisition, the correction should be applied as the first 
step in data reduction and analysis. 

As suggested by ^Bi retta fc Kozhurina-Plataij (|2005h . 
warm pixels in the detectors were successfully used to mea- 
sure device properties including the charge trap density and 
release times. The ubiquitous distribution of warm pixels 
allows the locations of charge traps to be continuously mea- 
sured across the CCD, using in-orbit science data. Exceed- 
ing earlier expectations, however, the range of fiux levels in 
warm pixels can also be used to obtain the 3D height of the 
traps, and the rate at which electrons fill up the volume of 
a pixel's potential well. All of these parameters are needed 
to correct CTI trailing, and it is remarkable that they can 
all be measured using a separate consequence of the radi- 
ation damage. To improve on this with future large sur- 
veys will probably require laboratory-based measurements 
of the CCD and charge trap characteristics before launch 
(and possibly after retrieval), or the integration of on-board 
electronics to allow in-fiight tests such as pocket pumping. 

We find two significant species of charge traps that af- 
fect the astrometry and morphology of astronomical objects, 
with characteristic release times r = 10.4 ± 3.2 pixels and 
T = 0.88 ± 0.2 pixels. These accumulated in a ratio of ap- 
proximately 3.0:1 and, by the end of HST cycle 13, there was 
slightly more than one in every other pixel, uniformly dis- 
tributed throughout the CCD. An electron beginning 2048 
pixels away from the readout register therefore encounters 
a large number of traps during readout. The (random) loca- 
tions of traps in our model are inevitably different to those 
in the real CCD. This adds shot noise during correction 
that is particularly noticeable in the structure of the image 
background. We have eliminated half of the shot noise by 
increasing the density of model traps density and allowing 
them to hold non-integer amounts of charge. 

We did not attempt to measure th e density of traps wit h 
very long release times. As shown bv [Rhodes et al\ l|2009h . 
these primarily affect object photometry. They can be mea- 
sured using repeated images of standard stellar fields, such 
as the STScI internal CTE calibration programme. Alter- 
natively, once the uniform distribution of other charge trap 
species has been established, they could also be measured 
from EPER or FPR tests. If their density is measured after 
servicing mission 4, our software can easily handle additional 
species of charge traps, and thus improve the simultaneous 



correction of object photometry, astrometry and morphol- 
ogy- 

We measure a functional form describing the rate at 
which electrons fill up the volume of a pixel's potential well. 
Below a "notch" depth of 96.5 ± 2.0 e~, they occupy neg- 
ligible volume; additional electrons overflow and the cloud 
expands with a power law of index a = 0.576 ± 0.013. This 
volume determines the number of traps encountered by the 
cloud, as it is moved through the CCD substrate towards 
the readout register. Note that, in order to cope with the 
huge data volume, we developed the model stepwise: fix- 
ing the first measured parameters before fitting the later 
ones. Since errors on some measurements were not carried 
through, these later errors are likely to be underestimated. 

Our pixel filling model neatly illustrates why CTI trail- 
ing is a non-linear function of object flux, and why it can 
be reduced by pre- or post-flashing an image to increase the 
background level. During readout, the large cloud of elec- 
trons within a bright source has to navigate more charge 
traps that the small cloud within a faint source - but the 
ratio of exposed traps to electrons falls as Uq/ue ~ n"~^ . 
Since q < 1, the fractional amount of charge that gets trailed 
behind a faint source is greater than that behind a bright 
source. In addition, a fat zero reduces the total amount of 
trailing behind all sources by moving all pixels to the right 
in figure 13] For faint sources of fixed fiux, the total number 
of exposed traps is proportional to Auq/Aue, which decreases 
as ~ {b~dy-^. 

The two biggest simplifications in our model are reduc- 
ing the three-phase clocking cycle to a single shift, and not 
explicitly treating the continuous injection of charge from 
warm pixels during readout. The latter would be particularly 
important if this method were required for measurements of 
charge traps with long release times. Improvements in these 
areas would be interesting directions for future studies. 
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